clear all 
clc
y=[63.1 63.9 65.6 66.8;65.1 66.4 67.8 69.0;67.2 71.0 71.9 73.5];
p=anova2(y)
